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ABSTRACT 

We present results from the first application of the Grid of i?ed Supergiant and Asymptotic Giant 
Branch Models' (GRAMS) model grid to the entire evolved stellar population of the Large Magellanic 
Cloud (LMC). GRAMS is a pre-computed grid of 80 843 radiative transfer (RT) models of evolved 
stars and circumstellar dust shells composed of either silicate or carbonaceous dust. We fit GRAMS 
models to ~30 000 Asymptotic Giant Branch (AGB) and Red Supergiant (RSG) stars in the LMC, 
using 12 bands of photometry from the optical to the mid-infrared. Our published dataset consists of 
thousands of evolved stars with individually determined evolutionary parameters such as luminosity 
and mass-loss rate. The GRAMS grid has a greater than 80% accuracy rate discriminating between 
Oxygen- and Carbon-rich chemistry. The global dust injection rate to the interstellar medium (ISM) 
of the LMC from RSGs and AGB stars is on the order of 1.5 x 10~ 5 M Q yr _1 , equivalent to a total 
mass injection rate (including the gas) into the ISM of ~ 5 x 10 -3 M yr _1 . Carbon stars inject two 
and a half times as much dust into the ISM as do O-rich AGB stars, but the same amount of mass. 
We determine a bolometric correction factor for C-rich AGB stars in the K s band as a function of J 
- K s color, BC Ks = -0.40(J - K s ) 2 + 1.83(J - K s ) + 1.29. We determine several IR color proxies 
for the dust mass-loss rate (M^) from C-rich AGB stars, such as log Md = ^ K _"j 8 18 Q^ 3 37 — 5.93. We 
find that a larger fraction of AGB stars exhibiting the 'long-secondary period' phenomenon are O-rich 
than stars dominated by radial pulsations, and AGB stars without detectable mass-loss do not appear 
on either the first-overtone or fundamental-mode pulsation sequences. 



1. INTRODUCTION 

At the end of their lives, stars of approximately solar- 
mass (0.8-8 Mq) ascend the Asymptotic Giant Branch 
(AGB), the final phase of nuclear burning in the lives of 
these stars. The AGB is one of the brightest populations 
in the infrared (IR) sky, contributing up to ~20-30% 
of the IR light for the Small Magellanic Cloud (SMC) 
(Boyer et al. 2011). Essentially all AGB stars are vari 



ables ( Vijh et al 



hundreds ot days. 



2009), with periods on the order of 
his variability is caused by hydrody- 



namic pulsations traveling through the extended atmo 
sphere of the star. The dramatic changes in stellar ra- 
dius these shocks produce c ause brightness flu ctuations 
on the scale of ^2 mag (e.g. Wood et al.||1999|), on time 
scales of hundre ds of days (Eraser et al. 2008 White- 
'Riebel et al.|[201U[). On much longer 



lock et al.||2008[ 
time scales (~10 yr), runaway thermonuclear reactions 
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in the helium burning shell of an AGB star, called 'ther 
mal pulses, ' can cause much more dram a tic brightness 



variations (|Schwarzschild fc Harm 1965 Vassiliadis & 
Wood|1993p . Thermal pulses cause readjustments in the 
glo bal st ruct ure of the AGB star, carrying elements syn- 
thesized in the nuclear burning regions to the surface. 
This process, known as the 'Third Dredge Up,' is respon- 
sible for the formation of carbon-rich AGB stars ( |Iben| 



19831, as well as b ringing more exotic elements su ch as 
I'c to the surface (Uttenthaler & Lebzelter 2010), and 
its precise details re main a topic of current research (e.g. 
|Karakas et al.|20l"o" ) . Towards the end of their e volution, 
AGB stars exhibit extensive rates of mass-loss ( |Wachterj 
|et al.|2002[ ) , driven by the aforementioned hydrodynamic 



pulsati ons and radiation pressure on the resultant dus t 



grains ( |Winters et al.||2000| |Mattsson fc Hofner||2011 ) 
This mass-loss, enhanced by elements produced in the 
nuclear-burning regions of the star and dredged to the 
surface, makes AGB stars one of the primary sources for 
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a-elements in the universe, and an important contribu- 
tor to dust in galaxies. The precise degree to which AGB 
stars contribute dust to the Interstellar Medium (ISM) is 
a matter of some dispute however. Specifically, to what 
extent do AGB stars or supernovae (SNe) dominate this 
process? 

The AGB contribution to a galaxy's total dust budget 
can be computed if the entire population of mass-losing 
AGB stars in the galaxy is identified, and the rate of 
dust production by each star in this sample is known. A 
comprehensive study of Galactic AGB stars is hampered 
by extinction in the plane of the Milky Way. The Large 
Magellanic Cloud (LMC) offers the ideal combination of 
relative proximity and low line-of-sight extinction, allow- 
ing for detailed studies of galaxy-wide evolved star pop- 
ulations. The LMC's high galactic latitude minimizes 
both foreground contamination by Milky Way stars and 
reddening due to intervening dust. In addition, the dis- 
tance to the LMC (~50 kpc) is well d e termi ned (e.g. 



Table 1 

De-reddening 
Coefficients 



Band 


A x (mag) 


U 


0.5900 


B 


0.5315 


V 


0.4590 


i 


0.2708 


.1 


0.1125 


H 


0.0652 


K s 


0.0372 



Ngeow & Kanbur 20081 



(2008) for dis- 



See Schaefer 

cussion and meta-analysis of this measurement. This is 
close enough that individual stars can be resolved, yet far 
enough away that the 3-D structure of the LMC can be 
neglected, accurate distances assumed for all stars, and 
therefore intrinsic brightnesses determined. 

Recent large-scale photometric surveys of the LMC 
such as the Magellanic C louds Photometric Survey 
(MCPS; | Zaritsky et aL][2004|), the Tw o Micron All-Sky 

and the Survey- 



Note. — De- 
reddening coefficients 
used in this work for 
non-SAGE photome- 
try. SAGE pho- 
tometry was not de- 
reddened, as interstel- 
lar reddening is neg- 
ligible at those wave- 
lengths. 

rich and carbonaceous dust respectively. In this paper we 
apply the GRAMS to the whole population of evolved 
stars in the LMC in order to derive more precise mea- 
surements of individual sources and the population wide 
return of mass to the galaxy. 

The remainder of this paper is organized as follows: 
Section [2] of this p aper discusses the sources of our ob 



Survey (21V1ASS; |Skrutskie et al.||2006 
inq the Agents of a Galaxy's involution survey (SAGE; 
Meixner et al. 2006), have allowed the construction of 
catalogs of the entire AGB population of that galaxy, 
with mu lti-band photometr y for tens of thousands of 
sources dBlum et al.|[2 006|). The SAGE-Sp ec follow 



se rvat ional datajjj pO) and our fitting procedure (§ 2.2 
& |2.3[ ) . Section (^displays some specific fits representa- 
tive of our sample. Section [4] describes the major results 
of our fitting, including the O-rich/C-rich determination 
(§ |4.1[ ), the revealed luminosity fun ction of the evolved 
stellar population of the LMC (jj]4.2|), the integrated dust 



mass return to the ISM (§ 4.3 &FTA), and the develop 



on to the SAGE survey lKemper et al.| ( |2010[ ) has pro- ment of simple observationaTproxies for dust mass-loss 

_1_.___._1 _.._____..___! _l__._. , r , ___ ^^^^^^^^^^^^™^^^^T^^^^^^^^ . . J 1 _ _ T 1\ ff ~v f\ 4~ s~i fK A ^ 1 \ A / _-\ vat T i attt _-i -n nlllfll n t v\ _7I 1~\ r i f-^ -i-i 1^1 



duced spectral class ifications tor H7D so urces m the LMC 
( |Woods et al.|2011| The IRSF survey (iKato et al.|2007 l 
also examined the LMC in near-IR (,1HK S ) bands. 



rate (§ 4.5). We review our conclusions in section [5j 



The AKARI mission (Murakami et al. 2007 ) has also 
been used to survey a portion of the LMC (Ita et al. 
20081. While this survey covers a smaller area, the 
unique [11] and [15] /im bands of the AKARI satellite are 
a valuable contribution to the multi-wavelength coverage 
of this important gala xy. The ongoing V ista Magellanic 
Clouds Survey (VMC) |Cioni et al.|(l2011[ ) will push three 
magnitudes deeper than the AKARI survey, and has al- 
ready been used to develop a spatia lly resolved star for - 



mation history of part of the LMC (Rubele et al. 



As a first step towards the direct measure ment o 



2011) 



mass- 



los s rates for the entire LMC AGB sample, (Srinivasan et 
al.||2009l hereafter Paper I) computed infrared excesses 



2. DATA AND FITTING PROCEDURES 
2.1. Data 

The original SAGE survey w as conducted in tw o 
epochs, spaced ^3 months apart ( Meixner et al.| [2006). 
The observations from these epochs have been combined 
into a single mosaic photometry archive and catalog, 
which is deeper and has smaller photometric errors than 
the individual epochs. Our dataset consists of 33 718 
sources extracted from the SAGE Mosaic Photometry 
Archive. These sources have been matched to optical 
data fro m the Magellanic C louds Photometric Survey 
(MCPS) ( iZaritsky et al.|2004|), near-infrared ph otometry 
from the 2MASS survey ( fsKrutskie et al .||2006|), and the 



from SAGE data and used these to estimate the total 
dust injection rate. Boyer et al. (2012, in press) have 
performed a s imilar analysis for the Small Magellanic 
Cloud (SMC). patsuura et al] ( |2009| ) determined a sim- 
ple IR col or proxy for mass-loss us ing mass-loss rates 
derive d by Groenewegen et al.| (]2007b and|Gruendl et al. 



( 2008 ) through detailed modeling of individual sources 



and extrapolated a global gas and dust budget for the 
LMC. An alternative to detailed modeling is to compare 
the observed S EDs of sources to pr e-computed models. 
Wit h this aim, Sargent et al. ( 2011| hereafter Paper IV) 
and |Srinivasan et al.| ( |201l[ hereafter Paper V) presented 
the Grid of RSG and ATJEJModelS (GRAMS) for oxygen- 



variability information of the MACHO survey ( Alcock et 
al. 1999[ ), allowing us to construct 12 band spectral en- 
ergy distributions (SEDs) from the U band to 24 /zm for 
most of our sources. A 2" matching radius was used for 
all catalog joins. Details of the matching to the MCPS 
catalog can be found in the SAGE Data Delivery Doc- 
umentR The join to the MACHO cat alog is discussed 
more thoroughly in [Riebel et aL] ( |2010[ ). All non-SAGE 
photometry has been de-reddened, with the de-reddened 
flux, F , related to the observed photometry, F obs by 



1 http://data.sp_.tzer. caltech.edu/popular/sage/ 

20090922_en\discretionaryi->i>t>_ianced/docu__ents/ 
SAGEDataProductsDescription_Sep09 . pdf 
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Fo = -Fobs x 10( 0AAx \ Our de- reddening coefficients, A\, 
are listed in Table [TJ 

When extracting sources from the S AGE database, our 
initial AGB classifications follow from jCioni et al. (20061 
and Blum et al. (2006). A star is classified as an oxygen- 
rich (U-richj or carbon-rich (C-rich) AGB candidate 
based on its location on the K s vs. J — K s color magni- 
tude diagram (CMD). The J — [3.6] color (or, in the ab- 
sence of a J-band detection, the [3.6] — [8.0] col or) is used 
to s elect extreme AGB candidates. See § 2.2 in |Riebel et| 
al. (2010) for explicit definitions of these color cuts. RSG 



candidates are sel ected based on the color-magnitude cri- 
teria presented in Boyer et al. ( |2011 ). Specifically, RSGs 
candidates are defined as being brighter than the Tip 
of the Red Giant Branch (TRGB), K s = 12, and be- 
tween the lines K s = -13.333(J - K s + 0.25) + 24.66 
and K s = -13.333(J - K s ) + 24.66. These CMD-based 
definitions are illustrated in Figure [T] 



Looking fo r Young Stellar Objects (YSOs), [Gruend l 
~" (|2CI08|) published a list of 13 'Extremely Red Ob 



et al 



jects' that they spectroscopically identified as carbon 
stars, based on the presence of SiC absorption. Twelve 
of these sources lie outside our defined CMD cuts and are 
therefore not included in our initial dataset, but we have 
manually added them to our list, using the photometry 
published in that work. All of these sources are classified 
a s extreme A GB candidate s. 

|Paper IV] and |Paper V present the development of 
GRAMS, a grid of RT models of dusty evolved stars cal- 



culat ed using the 2-Djust RT code (Ueta & Meixner 
20031. The dust properties for thes e models were fixed 



by modeling two O-rich AGB stars ( Sargent et al.|[2010 



hereafter Paper II) and one C-rich star in the LMO 



(Srinivasan et al. 



2010 



hereafter Paper III). GRAMS 
consists of 68 600 O-rich and 12 243 C-rich models, span- 
ning a large parameter space of stellar photosphere and 
circumstellar dust shell properties. The output from 
GRAMS consists of spectra as well as synthetic pho- 
tometry for a large set of narrow- and broadband filters, 
including the ones used in our study. Using these pre- 
computed SEDs, we are able to find the best fit model to 
each source in our sample through a simple brute force 
search on an average desktop computer in only 4 hours. 
This approach may reduce the detailed accuracy of any 
particular model fit, but will compensate by allowing sta- 
tistically accurate trends and patterns to be determined 
for entire stellar populations in computationally reason- 
able periods of time. Figure [2] shows the entire GRAMS 
grid on a [8.0] vs. [3.6] - [8.0] CMD. The O-rich grid 
is shown in the left panel and the C-rich grid is shown 
in the right panel. The gray background is the entire 
extent of the grid, and the highlighted points (blue for 
O-rich, red for C-rich) emphasize the models that have 
been matched to sources in our dataset. 

2.2. Variability 

The overwhelm ing majority of AGB stars are variable 
( Vijh et al.|2009), with co mplicated, multi-periodic light 
curves ( |l'raser et al.|2008[ ) . Because our dataset combines 
multi-epoch observations from different surveys taken on 
uncorrelated dates, the optical, near-IR and mid-IR por- 
tions of the observed SEDs for our stars sample different 
points in the stars' light curves. Accurate phase correc- 
tion for such complicated light curves for a sample the 



size of ours is not practical, nor possible for the many 
sources which lack variability data from long-term mon- 
itoring surveys like MACHO. As a major goal of this 
project was to demonstrate large-scale radiative transfer 
modeling, and the fact that many of the sources without 
MACHO data are the extreme AGB candidates which 
dominate the mass-return to the ISM we seek to measure 



(Paper I), we opted to handle the sources' variability in 



a statistical manner rather than restrict our sample to 
those objects for which precise phase-correction could be 
determined. 

To account for variability, we took the SAGE mid-IR 
photometry as the baseline observations, and inflated the 
error bars on the flux measurement in the shorter wave- 
length bands (UBVIJHK) by adding an additional error 
term representative of the source's variability amplitude 
to the photometric errors. Beca use the MACHO su rvey 
used a non-standard filter set (Alcock et al. (19991), we 
first transformed the MACHO 



1999[ ), 

3- band photometry into 



a Johnson V magnitude using the prescription in 
et al.l (119991), 



ry 
Al 



lcock 



V = 24.114 + 1.00258b - 0.153(b - r) 

where b and r are the MACHO mean b- and r-band mag- 
nitudes, respectively. This transformation allowed us to 
easily cast the MACHO variability amplitude (originally 
in magnitude units) as a ratio of the source's flux at min- 
imum brightness to that at maximum brightness. We 
choose the MACHO b-band because it is typically larger 
in amplitude than the r band, and thus we err on the 
side of caution when modifying our error bars. 

We modeled the sources as being a constant flux source 
with a single period sine- wave signal imposed on it. We 
used the b-band mean MACHO magnitude as the aver- 
age brightness of the source, due to the seven-year base- 
line of this measurement. The amplitude of the variation 
imposed on this mean was the MACHO b-band ampli- 
tude assoc iated with the dom inant period of variation, 
taken from Fraser et al. ( 2008 ) . It can be shown that the 
RMS average of a sine function imposed on a constant 
flux Fy is 

(\-u 

Cvar — Fy 



1 



a 



where a is the ratio of the minimum to maximum flux, 
i.e. the MACHO amplitude (measured in magnitudes) 
converted to flux units. One-half this quantity was added 
in quadrature to the photometric errors of all non-Spitzer 
bands for each source with available MACHO variability 
information. The factor of one-half accounts for the fact 
that the variability amplitude of AGB stars is greater 
when measured in the optical (i.e. at the wavele ngths 
studied by MACH O) than in the near-IR (see, eg. Reid 
fc Goldston||2002| ) . For those sources in the SAGE cat- 
alog without MACHO observations, a "canonical" vari- 
ability amplitude was constructed from the sources which 
did. All sources were classified as O-rich, C-rich, or ex- 
treme AGB c andi dates, using the photometric criteria 
described in § |2.1| The median amplitude of each class 
was then assigned to all stars of that class without MA- 
CHO data, and we inflated the error bars of the non- 
SAGE photometry in the same manner discussed above. 
That is, the median variability error term of all MACHO- 
detected sources classed as O-rich AGB candidates was 
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j - k s j - K s 



Figure 1. K s vs. J — K s CMD illustrating the two classification schemes used in this paper. Color-magnitude cuts (see § 2,2, |Riebel et| 
|al.|2010| for details) are first used to select AGB and RSG candidates (left panel). The observed SEDs of these candidates are then fit witn 
UKAMS models, which point to either oxygen-rich or carbonaceous dust chemistry (r ight panel). As the GRAMS chemical classification 
does not distinguish between O-rich AGB stars and RSGs, we use the CMD cuts of |Boyer et al. | ( |2011[ ) (thick lines in both panels) to 
identify RSG candidates. Sources with O-rich dust (blue dots) that fall within the CMD cuts are classified as RSGs. 




GRAMS O-rich Grid 



[3.6] - [8.0] 



2 4 6 
GRAMS C-rich Grid 



Figure 2. [8.0] vs. [3.6] - [8.0] CMD of the entire GRAMS model 
grid. Left panel: Every model of the O-rich GRAMS grid (light 
gray) with models selected as a best-fit to one of our sources high- 
lighted in blue. Right panel: Same plot for the C-rich GRAMS 
grid, with used models highlighted in red. The two grids were de- 
signed to span a larger area of parameter space than real stars are 
expected to cover. This plots shows that best-fit models are clus- 
tered in both grids and that the grids cover more CMD space than 
necessary. 



used as the "canonical" variability error term for all O- 
rich AGB candidates without MACHO detections, and 
similarly for C-rich and extreme AGB candidates. RSGs 
are not typically as variable as AGB stars, and we do not 
inflate the errors of RSGs candidates. 

2.3. Fitting Procedure 

The best-fit GRAMS model for each of our 33 718 can- 
didate AGB and RSG stars was found using a brute-force 
minimum % 2 search. Each source was compared to all of 
the -68 000 O-rich GRAMS models, and the -12 000 C- 
rich GRAMS models. The best-fit model was defined to 
be the model with the smallest value of the quantity 



N ^ 



where / bs ; and f mi are the observed and model flux in 
the i th band, respectively. N is the number of bands 
for which a source has valid photometry, and er$ is the 
quadrature sum of the photometric extr actio n error and 
the variability error term described in § |2.2| This quan- 
tity is thus properly a \ 2 per data point. For the red- 
dest sources, defined as those having a J/I flux ratio of 
10 or higher (equivalent to (I — J) > 1.4), we obtained 
better results by neglecting their optical (UBVI) pho- 
tometry entirely, which represents a negligible fraction 
of the energy of their SEDs. The J/I flux ratio was se- 
lected because this bridges the gap between the MCPS 
and 2MASS surveys, and is thus sensitive to both sources 
which are intrinsically very red and variable sources ob- 
served by the two surveys at very different points in their 
light curve. A star is classified as O- or C-rich based 
on whether the best fitting model is from the O- or C- 
rich grid. RSGs are not a separate classification within 
GRAMS, but the O-rich grid is designed to cover stars of 
higher luminosities than the classical AGB limit, includ- 
ing RSGs. When RSGs are discussed in this work, their 
cl assification as such is purely based on the CMD criteria 
of |Boyer et"aT1 ( |2011[ ) (illustrated in Figure [I]). Figure [3] 
shows the distributions of the x 2 values for model fits of 
both types. The median of the O-rich chi-squared distri- 
bution is much lower than that of the C-rich one. The 
differences in these distributions are taken into account 
when we address the reliability of classifications made 
by comparing the chi-squared values of the C-ric h and 
O-rich best fit models to a given source (see 
Figure [3] shows two typical fits to our dataT 



4.1}. 

These 



obsi 



fmi ) 



fits are 'typical' in that they were specifically selected to 
have values of \ 2 closest to the median value for their 
class (O- or C-rich). The left column shows the SED of 
the median O-rich fit above its location in the K s vs. J 
— K s CMD. The right column shows the same plots for 
the prototypical C-rich fit. Both of these sources are 
centrally located in their respective populations in the 
IR CMD (see Figure [l}. The effects of inflating the non- 
SAGE photometric error bars to account for variability 
can clearly be seen in the U- and B-band fluxes in the 
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Table 2 

Dataset Populations 



Figure 3. Distribution of the \ 2 P er data point for the best-fit 
model to each source in both the O-rich (blue) and C-rich (red) 
populations in our sample. Both distributions are sharply peaked, 
with long tails extending to large values. The C-rich sample in 
general has larger values of \ 2 f° r the best-fit models. 



C-rich SED, right column. 

After this initial fitting, we visually inspected the SEDs 
of several hundred sources. The distribution of x 2 for 
both the O-rich and the C-rich samples were found to 
be sharply peaked (Figure [3| , and we selected the 5% of 
the sample with the largest value of x 2 f° r inspection. 
For these stars we visually inspected the fit derived us- 
ing all valid photometry and the fit to the IR (2MASS 
and SAGE) data only. If the IR fit was significantly bet- 
ter (both quantitatively as defined by the value of x 2 
and qualitatively as defined by eye) the source was man- 
ually flagged to be fit using only the IR bands. Based 
on both manual identification and color criteria (J/I flux 
ratio), 6089 of our sources were fit using IR data only. 
We flagged 227 sources as invalid fits. 

The discarded sources included those with too few valid 
data points (< 4 bands), obvious foreground sources, and 
sources with SEDs not consistent with an AGB star. Fig- 
ure [5] displays some example sources with fits that were 
manually adjusted or rejected. At top left we show a 
source with an anomalously dim /-band flux. There were 
202 such sources in our dataset. Because the fit based 
only on the IR data points (red curve) matches the over- 
all SED much better than that based on all the photom- 
etry including the /-band (green curve), these sources 
were retained in our dataset, using the fit based only 
on the IR data. At the top right of Figure [5j we have 
a source which is obviously a foreground source, mis- 
identified in our catalog as belonging to the LMC. The 
best fitting model (blue curve) has fluxes ~10 times dim- 
mer than the observed object, yet a luminosity of a mil- 
lion L . If this source were truly in the LMC, it would 
be unphysically luminous. At the bottom right, we see a 
source with poor data quality. Such bright 2MASS pho- 
tometry compared to such dim [3.6] and [4.5] /mi pho- 
tometry is probably the result of a mismatch in the join 
between the 2MASS and SAGE catalogs. Such sources 
were removed from our sample. Finally, at the bottom 
left of Figure [5] we see a source with an SED not typical 
of an AGB star or RSG. The anomalous optical pho- 
tometry indicates that this source is either a mismatch 
between catalogs, or not an evolved star. 





Class (CMD) 


Class (GRAMS) 


0-rich a 


19 566 


26 210 


C-rich AGB 


6709 


7281 


Extreme AGB 


1340 


N/A 


Red Supergiant 


5876 


N/A 



Note. — Numbers of sources classified as each type 
(O-rich, C-rich, Extreme AGB, and RSG) using both 
color-magnitude cuts (see |Riebel et al.|2010| for details) 
and using the GRAMS model grid. CRAMS classifica- 
tions arc based on the class of the single best-fitting 
model, as defined by the smallest value of x 2 ■ 
a The GRAMS O-rich classification includes RSGs, as 
well as AGB stars 



We are left with a dataset of 33 491 sources matched 
to GRAMS models from which we obtain astrophysical 
parameters such as dust mass- loss rate (Md), luminosity, 
and effective temperature (T c ff). Table ^ summarizes 
these numbers. Our entire 33491 source catalog is avail- 
able online as a machine readable table. Table [3] presents 
a small excerpt as a guide to the structure ofthe elec- 
tronic table. For each source, we list the SAGE mosaic 
photometry ID, the classi fication assigned to the source 
using the CMD criteria of Cioni et al.| ( |2006| ), the classi- 
fication assigned by GRAMS, the \ z °f the best-fitting 
model, the M^, the luminosity, the T e g and the optical 
depth (r) of the circumstellar dust shell (determined at 
10 /zm if a star is classified O-rich, 11.3 /Ltm if C-rich) 
of that model, along with the associated uncertainties. 
MACHO variability data (amplitude and period) are in- 
cluded for sources for which it is available, along with the 
UBVIJHK, [3.6], [4.5], [5.8], [8.0], and [24] /tm fluxes (and 
uncertainties) for every source. The UBVIJHK photom- 
etry listed in the online table are not dereddened. Photo- 
met ric u ncertainties have not been inflated as described 
in § [2~2l 

The single best-fit model to a given source does not 
have uncertainty in its parameters such as luminosity, 
Md, and T c g (strictly speaking, M^ is not a parame- 
ter of the model, but it is an immediate consequence of 
optical depth, which is). In order to define an uncer- 
tainty for these parameters, we use the median absolute 
deviation (MAD) for each parameter from the 100 best 
fitting models in the same grid (O-rich or C-rich) as the 
uncertainty for that parameter. That is, for any given 
model parameter X (M^, luminosity, etc.), the quantity 
ax = Median(|Ai — Median(Ai)|), where i ranges over 
the 100 best-fitting models, is calculated and defined as 
the uncertainty in that parameter. The MAD is a more 
robust estimate for the spread in parameter values than 
the standard deviation or a similar statistic, and more 
appropriate for our purposes since it is highly unlikely 
that all of the output parameters from the best-fitting 
100 models to a given source will be normally distributed. 
Note that in the case of normally distributed errors, the 
MAD is smaller than the standard deviation, specifically 
stddev 1.5(MAD). For some ofthe more coarsely sam- 
pled parameters, such as T e ff, all 100 best-fitting models 
may have identical values of a given parameter. In these 
cases, the MAD will be 0, and in Table [3] we have set the 
uncertainty to the unambiguous placeholder value —99. 
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Figure 4. Two example fits from the GRAMS grid. These particular sources were selected because the quality of their fit (x 2 ) more 
closely matches the median value of x 2 f° r their respective grid than any other source. That is, they are the most 'typical' fits in our 
sample. Left panel: An O-rich source. Right Panel: A C-rich source. The bottom row shows the location of each source on the K B vs. J 
- K s CMD. 




Figure 5. Four sources selected for visual inspection of the fit. Top left: A source with an anomalously dim I-band flux. 202 sources had 
similar SEDs. The actual photometry is plotted in black, the best-fit model based on all valid photometry is in green, and the fit to only 
the IR data (2MASS and Spitzer) is in red. Notice the IR-based fit does a better job of matching all the available photometry because it 
is not skewed by the I-band. This source was kept in our sample, using the fit to only the IR data. Top Right: The unphysical brightness 
of this source indicates it is not at the distance of the LMC. This source was removed from our dataset. Bottom Left: A source with an 
SED inconsistent with an AGB star. The monotonic increase of the flux of this star as one moves to longer wavelengths is not typical of 
AGB stars, which peak in the near-IR. This source was removed from our dataset. Bottom Right: A source with bad data quality. The 
bizarre SED of this star is most likely due to a mismatch between the optical and IR surveys used in this study. This source was removed 
from our dataset. 
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Table 3 

Source List 



SAGE ID a 


Class b 


Class c 


x 2 






Lbol 


0"L 


Amplitude 


Period . . . d 




(CMD) 


(GRAMS) 




(M yr- 1 ) 


(Mg yr" 1 ) 


(Lq) 


(L©) 


(mag) 


(days) 


J050115. 85-692040.4 


o 


o 


5.56 


5.71E-11 


1.19E-11 


4149 


299 


0.13 


783.70 . . . 


J050311.47-691412.3 


o 


c 


17.36 


4.22E-11 


1.69E-11 


10000 


879 


0.08 


110.61 


J051227.54-701730.7 


X 


c 


42.53 


4.22E-10 


5.69E-11 


5128 


384 


1.74 


356.63 . . . 


J060647.79-664812.5 


s 


o 


31.52 


1.33E-08 


9.16E-09 


1000000 


-99 


0.00 


0.00 


Note. — Source list of datasct 


used in this s 


tudy. The full dataset is 


available in 


the online materials 


of this paper. 


This excerpt is 



included only as a guide to format. 

a In the online table, all SAGE-IDs are prefaced wi t h 'SST TSA OEMA' 

b CMD Classifications are described in |Cioni et al.| < |2006) and |Riebel et al.|j201o) 
c O-rich or C-rich classification based on best fitting GR AMS model 

The online data also contain the optical depth, T e ff, UBVIJHiC s , [3.6], [4.5], [5.8], [8.0], [24] fluxes, and their associated uncertainties 
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-10 

Figure 6. log^cr^/A/^) vs. log(Af,j) for our sample. O-rich 
sources in blue, C-rich in red. O-rich sources to the left of the 
gr een line at log(M [ j) = —11.3 arc classified as low-M^ sources 
(§ 4.3 1 and are consistent with = 0. All of these low-M^ 



sources have cr^/M^ > 1 while only 6% of sources with higher 
values of do. 

The number 100 was settled on after experimentation 
with various values from 10-1000. The number 100 is on 
the order of 1% of the GRAMS grid for both C-rich and 
O-rich models, and ensures that the majority of sources 
in our sample will have well defined uncertainties (not 
—99) for most parameters. Smaller subsets do not sam- 
ple a large enough region of parameter space, and often 
have all models in the subsample with identical values for 
many parameters. Larger subsets sample poorly fitting 
models and can cause an unacceptable number of models 
to have ax /X > 1. Using 100 models to define the uncer- 
tainty is a compromise between these two positions, with 
emphasis placed on the uncertainty in M^. The number 
of sources with valid errors for each parameter are listed 
in Table [4] We see that most sources have well-defined 
uncertainties for and luminosity. The uncertainty in 
T e ff is not well-determined from our model grid . Th is is 
an exp ected result. As discussed in |Paper IV] and Pa 



per V[ the model photospheres used in generating the 



GRAMS grid do not offer a very fine sampling of effec- 
tive temperatures. Figure [6] shows the logarithm of the 

ratio of the uncertainty in Md to the itself, log(-^-) 

vs. log(Aid)- O-rich sources with \og(M d ) < —11.3 (to 
the left of the green line in Figure pj) are not thought 
t o rep resent actual physical circumstellar dust shells (see 
§ 4.3 for discussion). Significantly, this value of M^ is also 
approximately where the S/N equals unity. Only 6% of 
sources with log(Md) > —11.3 have a signal-to-noise less 
than 1, while 100% of the low-M^ sources do. 



3. RESULTS: SED FITS AND TYPES OF SOURCES 

We fit GRAMS models to the SEDs of the O-rich stars 
(SSTISAGE1C J052206.92-715017.6 and HV 5715), as 
well as the C-rich star OGLE L MC LPV 285 79. These 
sources were modeled in detail in |Paper TI| and [Paper III] 
in order to establish the dust properties for the grid. Our 
results, listed in Table [5] agree to within uncertainties 
on almost every model parameter for all three sources. 
The exception, optical depth, is due to coarseness in the 



model grid coverage. The fact that our population scale 
fitting procedures can match the output of detailed mod- 
eling so well gives us a great amount of confidence in the 
overall accuracy of our results. 

Figures [7] and [8] show the SEDs and model fits to 8 
sources selected to illustrate the range of stellar parame- 
ters spanned by our sample. Specifically, we show sources 
with the maximum and minimum values of both M^ and 
luminosity for each class of source, C-rich and O-rich. In 
both figures, the right hand column shows an IR CMD 
highlighting the specific source shown in the left column. 

The top source in figure[7]was chosen because it has the 
largest value for M^ in our C-rich sample. Its SED peak 
is shifted far into the IR, peaking redward of 8 /im. It is 
a bright source, but only half as bright as the brightest 
C-rich sources in our sample (se e bel ow) . This is consis- 
tent with Fig 16 (discussed in § |4.3[ ), where there is no 
clear relationship between and luminosity seen in our 
sample for either C-rich or O-rich sources. This source is 
highlighted in the [8.0] vs. [3.6] - [8.0] CMD in the right 
panel. This Spitzer CMD is used because this source is 
so red that it is lacking a valid J or K s detection, and 
therefore cannot be shown on the J — K$ vs. K s CMD. 
As one would expect for such a red source with a high 
Md, this source is class ified as an extreme AGB star us- 

As discussed 



ing the color criteria of Blum et al 



in S 14.31 these extreme AGB stars dominate the mass re 



(2006). 



turn to the ISM in the LMC, accounting for 74% of the 
total amount of mass lost, but only 4% of the sample by 
number. 

The source shown in the second row of Figure [7J has 
the lowest M^ of our sample. It is one of the few stars in 
our sample classified as C-rich by GRAMS that qualify 
as a l ow-M(j star, defined as having \og(Md) < —11.3 
(§0). The uncertainty in its M d (±1.3x 10" U M Q yr^ 1 ) 
makes this source consistent with a bare photosphere and 
zero Md- The SED of this source peaks in the H band, 
consistent with a much hotter object (^2500 K) such as 
an actual visible star not surrounded by dust. 

The third and fourth rows of Figure [7] show sources 
with the maximum and minimum values (respectively) 
of luminosity for the C-rich grid. Both of these sources 
have similar near-IR colors and similar low values of M f /. 

Fig [8] has the same format as Figure [7] but all the 
sources here are classified as O-rich. Again, the source's 
location on an IR CMD is shown in the right column. 

The top source in Figure [8] has the highest M^ of any 
O-rich source in our sample. Its J — A" s color is extremely 
red because the SED peak has been shifted all the way 
to ^24 /im by thick circumstellar dust. We see a sharp 
contrast with the low-M^ source (essentially a bare pho- 
tosphere) seen in the second row. For the low-M^ source, 
the SED peaks in the H band. 

The bottom two rows of Figure M are the highest and 
lowest luminosity sources in our O-rich (which includes 
models of RSGs) sample. The lowest luminosity source 
shown at the bottom has a very low M^ as well, but the 
brightest source (1 x 10 6 L Q ) shows quite a high rate of 
luminosity driven mass-loss. 

4. DISCUSSION 
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Table 4 

Sources with Valid Uncertainties 



Class (GRAMS) 


Luminosity 


M d 


T e ff 


T 


O-rich 


25 579 (97.6%) 


26 210 (100%) 


7710 (70.6%) 


25 626 (97.8%) 


C-rich 


6963 (95.6%) 


7281 (100%) 


6996 (96.1%) 


6815 (93.6%) 



Note. — Number of sources with valid uncertainties for each parameter obtained from the 
model grid. Uncertainties in model parameters arc defined as the median absolute deviation 
of that parameter for the 100 best-fitting models to a given source. If all 100 models have 
the same value for a given parameter, the median absolute deviation is mathematically 
and a placeholder value of —99 is listed. 



Table 5 

Comparison of fitting results for sources also in |Paper II| and |Paper iH] 



ID b GRAMS Class M d (xlO" 9 M yr" 1 ) Luminosity (L ) T eff (K) 



SAGE ID a 



J051811.08-672648.5 d 


82 


O 


1.5 ± 


0.4 


33 694 ± 


6000 


3500 ± 400 


0.0256 


(HV 5715) 







2.3(1.1 - 


-4.1) 


36 000 ± 


4000 


3500 ± 100 


0.012 ±0.003 


J052206. 92-715017.7° 


96 





2.1 ± 


0.4 


4820 ± 


670 


3700 ± 200 


0.1024 


(SSTSAGE052206) 







2.0(1.1 - 


-3.1) 


5100 ± 


500 


3700 ± 100 


0.095(0.07-0.13) 


J051306.40-690946.3 f 


66 


c 


2.1 ± 


0.4 


7080 ± 


700 


3100 ± 200 


0.4 


(LPV 28579) 




c 


2.5(2.5 - 


- 2.9) 


6580(6150 


- 7010) 


3000 


0.27(0.25 - 0.275) 



Note. — Comparison of the stellar parameters for the sources examined in |Paper II| and Paper IIlJ Each source is listed twice, first 
with the parameter values found in this work, and second with the values found in the previous study in which it was modeled more 
precisely. If no uncertainty was found for a parameter (sec § |2.3| l , none is listed. Where uncertainties are not symmetric about the best 
fit value, the range is given in parentheses. 

a Source IDs from this work are all prefaced with 'SSTISAGEMA' 
b ID from the SAGE-Spec survey 

c The quoted optical de pth is tha t at 10.0 fim for sources in |Paper II| and at 11.3 fim for the : 

d Source HV 5715 from |PaperIl| 

c Source SSTSAGE05220 6 from|Paper II] 
f Source LPV 28579 from |Paper IIl| 

4.1. O-Rich/C-Rich Discrimination 



from [paper IlT| 



star would reach the Chandrase khar limit, based on the 



Cioni et al. (2006) proposed a photometric means of 
discriminating between O- and C-rich AGB stars in the 
K s vs. J - K s CMD. The GRAMS grid shows good 
agreement with these color-magnitude cuts, giving us 
confidence in the grid's ability to reproduce previously 
established population-scale classification schema . Ta- 
blc[2~l summarizes our agreement with the cuts from |Cioni| 



luminosity-core mass relation of |Paczyhski ( 1971 ) . How- 
ever, the most massive AGB stars can violate this limit 



due t o hot-bottom burning (Bloecker & Schoenberner 



1991 1, which also tends to prevent a star trom becom 



ing C-rich (Boothroyd et al. 



1993). 



et al.| (2006). Of the 26 210 sources we classify as O-rich, 
95% are classified as eit her O-rich AGBs using the cuts 
from |Cioni et al.| (|2006|) or as RSGs using the cuts from 



Boyer et al. J 2011 ). We also obtain a 83% agreement be- 
tween s ources classified as C-rich by the color-magnitude 
cuts of Cioni et al. (2006) and sources we match to C- 
rich models. The extreme AGB star classification does 
not exist as a separate category in the GRAMS model 
grid. One of the aims of this project was to provide 0- 
or C-rich classification for these sources, most of which 
are expected to be highly evolved and thus C-rich. Of 
the 1340 sources classified as extreme based on CMD 
definitions, 1299 sources (97%) are matched to a C-rich 
GRAMS model. Figure [9] shows the SEDs of two sources, 
both classified as extreme AGB stars by CMD criteria. 
The left column shows a star which is best fit by a C-rich 
model, primarily on the strength of its SAGE photom- 
etry, and the right column shows a star classified as an 
O-rich source, primarily because of its optical photome- 
try. 

We find 29 sources (not classified as RSGs via CMD 
cuts) in our sample with luminosities greater than the 
classical AGB limit, 54 000 L . This luminosity limit is 
based on the point at which the core mass of an AGB 



We classify a source as O- or C-rich based on the best- 
fitting GRAMS model. We assign a confidence in this 
classification by examining the ratio Xbest/Xaif wnere 
the two values of x 2 are the best fitting models from 
each GRAMS grid. Because the \ 2 distributions of the 
two classes of star peak at different values (Figure [3]), we 
use different criteria to define the confidence intervals for 
each class. For sources classified as O-rich, we consider 
them to be high confidence if xblxb < 0.1 and low con- 
fidence if xb/xb — 0-2 (medium confidence is interme- 
diate to these two) . For sources classified as C-rich, high 
confidence is considered to be xblxb < 0-4 and the clas- 
sification is considered low confidence if xblxb — 0<6. 
Tab le |6| gives the number of stars in each category. Fig- 
ure 1 1 Q| shows a histogram of these ratios for both cat- 
egories of sources. The vertical lines denote the points 
dividing high, medium and low confidence. The division 
points were selected to fall where the source density has 
fallen by a factor of approximately 2. 

We have tested our classification against the spectro- 
scopic ally v erified sources us e d by J Groenewegen et al"] 
(2009) and van Loon et al. (1999). All sources were 
matched to our catalog using a 2" matching radius. Of 
the 66 C-rich sources used in Groenewegen et al. (2009) 
that we find in our sample, we correctly classify 60 as C- 
rich. Of 27 C-rich sources in van Loon et al. (1999), we 
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Figure 7. C-rich sources showing the largest /smallest values of the and luminosity parameters found in our sample. The left column 
shows the SED fit for each source, and the right column shows its location in the K B vs. J — K B CMD (if the 2MASS photometry for a 
source is not available, the [8.0] vs. [3.8] — [8.0] CMD is shown instead). Top row: Highest value of M^. Second Row: Smallest values of 
M^. Third Row: Highest luminosity. Fourth Row: Lowest luminosity. 



classify 24 of them as C-rich as well. These numbers cor- 
respond to 91% and 89% agreement, respectively. The 
O-rich sources show s imilarly good agreement. We find 
40 of the sources that |Groenewegen et ah ( 2009 ) classify 
as O-rich in our sample, 38 of which we als o identify as 
O-r ich. O f the O-rich sources identified by |van Loon et| 



al 



(1 19991), 14 are found in our sample and 10 are identi 
rIed15yl~!RAMS as O-rich. 

In addition , we co mp ared our sample to the surve ys of 
|Blanco et al"1 ( [T980| and |Blanco fc McCarthy ] | ^ 1990 ] ) . We 
identify 177 carbon stars from Blanco et al.| (|l980 ) in our 
sample, correctly cl assifying 145 ( 82% ) of tnem. Of the 
96 O-rich stars from Blanco et al. ( 1980 ) identified in our 



dataset, GRAMS classified 92 (96%) correctly. |Blanco fc 



McCarthy ( 1990[ ) focused exclusively on C-rich stars, and 
we find 538 of their sources in our list, with 426 (79%) 
classified correctly by us. 
We also compare our sa mple to the point sou rces from 



the SAGE- Spec program ( Kemper et al.pOlO ) classified 
by|Woods et all (Ml). We find 87 sources from that 



project with valid GRAMS model fits in our own data, 
and excellent agreement between the two classifications. 
Of 17 SAGE-Spec RSGs, 100% are matched to GRAMS 
O-rich models, of which the RSG models are a subset. Of 
the other 70 SAGE-Spec sources, only 3 have conflicting 
GRAMS and SAGE-Spec classifications, a 96% success 
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Figure 8. O-rich sources showing the largest/smallest values of the and luminosity parameters found in our sample. The left column 
shows the SED fit for each source, and the right column shows its location in the K B vs. J — K B CMD (if the 2MASS photometry for a 
source is not available, the [8.0] vs. [3.8] — [8.0] CMD is shown instead). Top row: Highest value of M^. Second Row: Smallest values of 
M^. Third Row: Highest luminosity. Fourth Row: Lowest luminosity. 



rate. Two of the 3 sources with conflicting classifications 
are low-confidence as defined earlier in this section (Fig- 
ure 
AG 



12 top two panels). The remaining source (SSTIS- 
MA J053027.49-690358.3, Figure[l2j bottom panel) 



is technically a high confidence source, but examination 
of its SED shows that it is extremely well fit by both 
an O-rich (xo = 1-5) an d a C-rich model (xh = 5.9). 
The fits are essentially indistinguishable for this source. 
Combining all the above comparisons, we conclude that 
GRAMS has a greater than 80% accuracy rate when com- 
pared to spectroscopic classifications, correctly classify- 
ing 786 out of 948 spectroscopically classified stars to 
which our results were compared. The sources which are 



misclassificd by GRA MS lie very close to the O- and C- 
rich dividing line from Cioni et al. (20061 in the K s vs. 
J — K s CMD, squarely atop the locus of low-confidence 
classifications (green points, Figure 10). 
The 12 spect r oscop ically confirmed carbon stars from 



Gruendl et al. (2008) are all correctly identified as C- 
rich by GRAMS. We find that the current version the 
GRAMS grid does not contain models of sufficient optical 
depth to reproduce the 24 and 70 /im photometry given 
in that work, and we therefore use the luminosities and 
dust mass-loss rates derived by those authors. Beyond 
noting in passing that GRAMS did identify them as C- 
rich, these sources are only included in our dataset for 
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Figure 9. Two sources classified as 'extreme' AGB stars using the CMD criteria from |Cioni et al7| l |2006| ) ■ In both panels, the blue curve 
shows the best-fitting O-rich model, and the red curve shows the best-fitting C-rich model. Left Panel: A source classified as C-rich by 
the GRAMS model grid. 97% of the extreme sources are classified thus. Right Panel: A source classified as O-rich by the GRAMS model 
grid. The bottom row shows the location of the C-rich source in the K B vs. J — K B CMD, and the O-rich source in the [8.0] vs. [3.6] — 
[8.0] CMD due to the absence of JHK photometry. 



discussions of the global dust-mass injection into the ISM 
of the LMC. It should be noted that the mass- loss rates 



obtained by Gruendl et al. ( 2008 ) were derived using the 
sa me set of optical constants used by Groenewegen et 
al. (2007) which, as discussed in more detail in Paper V 



yield values of systematically 2-4 times higher than 
those produced by a GRAMS model of identical optical 
depth. 

Figure [Tl] presents a comparison betwee n the CMD- 
based classification scheme proposed by |Cioni et al.| 
(2006) and our current results. The top panel of the fig- 



ure presents the sources we classify as O-rich, with high 
confidence sources in blue, and low confidence sources in 
green. The bottom panel focuses on C-rich sources, using 
red for the high confidence sources and green again for 
lo w confidence sources . Both panels show the line defined 
by Cioni et al. (2006) which divides O-rich AGBs from 
C-rich AGBs in their schema. The GRAMS classification 
is in very good agreeme nt with the CMD-b ased classifi- 

is based on 



cation. As the line from Cioni et al 



(2006) 



only two bands of photometry, and GRAMS is based on 
12, we interpret this agreement as giving support to the 
CMD classification scheme. The fact that out of the 1337 
sources in this CMD with differing CMD/ GRAMS classi- 
fications, 1163 (87%) are considered low confidence clas- 
sifications lends support to our established confidence in- 
tervals. Where the G RAMS class i ficatio n disagrees with 
that derived from the Cioni et al. ( 2006 ) criteria, we give 
preference to the GRAMS classification. 

4.2. Luminosity Functions 

By covering the spectral region 0.2-1000 fim, the 
GRAMS grid provides accurate determinations of the 
bolometric luminosity of our sources. With this well- 
determined bolometric luminosity, we derive an expres- 
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Figure 10. Plot of the ratio Xbest/^ait f° r * ne sources classified 
as O- (top) and C-rich (bottom). Vertical lines denote the divisions 
between high, medium, and low confidence sources. 
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Table 6 

Confidence Intervals 



Class (GRAMS) 


Total 


High Confidence 


Medium Confidence 


Low Confidence 


O-rich AGB 


26 210 


16 609 (63.4%) 


4917 (18.8%) 


4684 (17.9%) 


C-rich AGB 


7281 


5213 (71.6%) 


1116 (15.3%) 


952 (13.1%) 



Note. — Number of classifications considered to be 'high,' 'medium,' or 'low' confidence. 
Confidence intervals are defined based on the ratio Xbest/x a lf where the two values of x 2 are 
those of the best fitting model, and the best fitting model from the other GRAMS grid (O- or 
C-rich). See text for details. 
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Ks) _ 

ble [7] We overlay the C-rich sample from |Gullieuszik "et] 



Figure 11. K B vs. J — K B CMD showing the sources classified 
as O-rich (top) or C-rich (bottom) sources by GRAMS with high 
(O-rich: blue, C-rich: red) or low (green, both panels) confi dence. 
Both panels feature the CMD- based division line proposed by |Cioni] 



(2006). For both classifications of star, the low confidence 
source are clustered near the CMD-based division. Of the 1337 
sources with different CMD and GRAMS classifications, 1163 are 
identified as low confidence. 



sion for the bolometric correction at K s (BCkJ as a 
function of J — K s color. We find a second-order poly- 
nomial fits our data well (Figure 13). Als o included 
in the figure is the quartic fit proposed by IWhitelock 
( 2006), and th e quadratic fit obtained by Ker- 
^ pOlO) . Our fit is of the form BC Kg = 
d2(J — K s ) 2 and is detailed in Ta- 



al.| (2012) as blue points atop our our dataset shown in 
reel. We note the remarkable consistency between these 
two samples. 

Figure 14 shows the luminosity function (LF) for the 



Table 7 

BC Kg vs. J — Kg 
Quadratic Fit 



Quantity 


Value 


an 


1.29 


ai 


1.83 


il-2 


-0.40 


N 


7281 


a 


0.06 



Note. — Fit pa- 
rameters for the best 
fit quadratic function 
to BC^ S as a func- 
tion of J — K B color. 
N is the number of 
sources used in cal- 
culating the fit, and 
a is the MAD of the 
residuals to the fit. 



RSGs are included in the O-rich sample. The maxi 



mum luminosity for an AGB star (^54 000 L Paczynski 
1971 1 is indicated with a vertical black line. The 'picket 



fence' effect seen at high luminosities is an artifact of 
increasingly coarse grid coverage at these high luminosi- 
ties (see also Figure [2]) , not a true characteristic of the 
RSG population luminosity function. The brightest O- 
rich source in our sample is a RSG far in excess of the 
maximum AGB luminosity limit, with a luminosity of 
1 000 000 L . The brightest C-rich AGB star in our sam- 
ple is 25 700 Lq , well below the maximum luminosity for 
a shell hydrogen/helium burning source. 

The O-rich distribution is concentrated at a lower lu- 
minosity than the C-rich distribution. C-rich sources 
are generally more evolved than O-rich sources, hav- 
ing had time to experience several thermal pulses and 
dredge up sufficient carbon to become C-rich. As lu- 



minosity increa ses with age along the AGB (Vassiliadis 
|fc Wood||1993[ though it should be noted that this is 
not a monotonic or smooth increase, but more of a 
'two steps forward, one step back' process), C-rich AGB 
stars should be brighter than O-rich AGB stars (see also 
However, the O-rich distribution does extend 
;o far higher luminosities than the C-rich popu 



Fig 117) 

a tail 



O-rich (black) and C-rich (red) populations of the LMC. 



lation reaches. The most massive AGB stars ignite the 
bottom of their convective hydrogen envelopes, and this 
'hot-bottom burning' ca n make them burn brighter than 
the classical AGB limit ( Bloecker fc Schoenberner|1991 ). 
Additionally, the AGB limit, obviously, does not apply 
to core-helium burning RSGs, which are also included 
amongst our O-rich grid. These two populations extend 
the bright tail of the O-rich LF, while the median of the 
C-rich LF remains brighter than that of the O-rich pop- 
ulation. 
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Figure 12. The 3 sources in our sample with different SAGE-Spcc 
and GRAMS classifications. Top/Middle Panels: These sources 
are classified by the GRAMS grid as low-confidence O-rich sources, 
a differing spectral classification is not surprising. Bottom P anel : 
This source is technically a high-confidence classification (§ |4.1| , 
but is fit excellently by both an O- and C-rich model. Again, a 
conflict between the spectrum and the RT model is not troubling. 



Figure [15] shows the integrated SED of the O- and C- 
rich populations of our sample individually, together with 
the combined total of our entire dataset. The figure was 
produced by simply summing the SED from each star's 
associated best-fit GRAMS model. This simulates the 
expected observed SED of the evolved stellar population 
of the entire LMC if it were unresolved. We note that the 
O-rich population (and its RSG sub-population) domi- 
nate the SED of the entire population at near-IR wave- 



Figure 13. Bolometric correction for C-rich AGB candidates at 
K B as a function of J — K B . The C-rich sample from this work 
is shown in red, and the best-fit second-order polynomial to the 
data is shown in green. For comparison, the C-rich sources from 
|Gullieuszik et al.| l |20f2| | are shown in as purple points. The fit 
is to all points in the figure. The dot t ed cu rve is the 4 th order 
polynomial derived in Whitclock et al. (2006). The violet dashed 
curve is the relation derived by [Kcrschbaum et al. (2 0f0| . 
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Figure 14. Luminosity functions for the stars in our sample. O- 
rich candidates are in black, C-rich candidates in red. The vertical 
line represents the classical AGB luminosity limit, ~54 000 Lq. 



lengths around 1 /im. The lesser number of C-rich stars 
means that they do not outshine the O-rich population 
(except by a small amount at mid-IR wavelengths), but 
do serve to "wash-out" the silicate feature at ~9 /im, 
making it much less prominent than if the O-rich stars 
were observed in isolation. 

4.3. Distribution of Dust Mass-Loss Rate Along the 

AGB 

One of the primary motivators for the study of AGB 
stars is their extensive mass loss. The broad picture of 



this p rocess is understood (see, e.g. Mattsson & Homer 



2011 1: hydrodynamic pulsations in the interiors of AGB 



stars develop, due to the steep density gradients in 
the extremely extended atmospheres of these stars, into 
shocks which eject the outer layers of the atmosphere. 
When this ejected material cools below the dust conden- 
sation temperature (^1000 K), dust forms. Due to its 
greater opacity, this dust is accelerated by the radiation 
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Figure 15. Plot of the integrated SED for the entire sample. Red 
curve: sum of all used GRAMS C-rich model SEDs, blue curve: 
sum of all used GRAMS O-rich SEDs, black curve: sum of red and 
blue curves. This figure constitutes a prediction for the integrated 
SED of the evolved stellar population of unresolved galaxies with 
metallicity and star formation history similar to that of the LMC. 
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Figure 16. M^ vs. bolometric luminosity for our sample. O-rich 
sources plotted in blue, C-rich sources in red. The left side of this 
figure is a scatter plot of all 30 000 sources in our sample. The 
discrete nature of the GRAMS grid leads to sources having best-fit 
models with identical values of M^ and luminosity, these sources 
are plotted atop one another in this figure. The left panel shows 
no clear trend of M^ with luminosity. The right hand panel of this 
figure is a histogram of M^ using the same color coding and scale 
as in the scatter plot. The O-rich sources are clearly bi-modal, 
and the minimum in this histogram (log(Af^) ~ —11.3) serves as 
the division between the low-M^ sources and the rest of the O-rich 
population. 



pressure from the star, and drags along the gas to which 
it is collisionally coupled. While it makes up only ~1% by 
mass of the material ejected from AGB stars, the dust 
takes on disproportionate scientific importance, due to 
its role as the driver of mass loss and because it is more 
easily observed in the infrared than the gas. The dust 
mass-loss rate is one of the stellar parameters spanned by 
the GRAMS model grid, and we have determined the 
for all ^30 000 sources in our sample. Throughout this 
work, we report only the dust mass-loss rate, not multi- 
plying by a gas to dust ratio, ip. Our reported mass-loss 
rates will thus be about 2 orders of magnitude smaller 
than values reported by studies which report tota l (gas 
and dust) mass loss (e.g Groenewegen et aL||2009| ). 
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Figure 17. Plot of [3.6] magnitude vs. variability period. C- 
rich candidates in red, O-rich candidates with log(Af ( j) < —11.3 in 
green, all other O-rich sources in blue. The 5 parallel AGB Period- 
luminosity sequences are visible (see also Figure 2, |Riebel et al.| 
1 2010 l). Low-M^ sources, thought to be the least evolved sources, 
are clustered on sequences 3, 4 and D. 



Figure 16 presents a plot of vs. bolometric luminos- 
ity, with a histogram of M<j presented along the right side 
y-axis. Luminosity and are presented on logarithmic 
axes, the histogram is linear. As both and luminosity 
are outputs of the GRAMS grid, the discrete nature of 
the grid is visible in this plot, particularly at the high lu- 
minosities. The fact that the grid is discrete also results 
in numerous "collisions," where multiple stars are fit with 
models with identical values of M^ and luminosity. These 
stars are plotted atop each other in the scatter plot, and 
are not visible as separate points. There is a slight ten- 
dency for Md to increase with increasing luminosity, but 
no quantifiable trend. For a given star, there is expected 
to be a trend of increasing M<2 with increasing luminos- 
ity during the course of its evolution. Since our sample 
consists of stars of various initial masses, various states 
of evolution, and various points along the thermal and 
hydrodynamical pulse cycles, this relation is washed out 
by the numerous complicating factors of AGB evolution. 
The Md histogram along the y-axis shows that C-rich 
stars are clustered at the high end of the mass-loss range, 
and the O-rich stars are extremely bimodal. We iden- 
tify a sub-population of O-rich sources based on this bi- 
modality, using M^ greater or less than log(Md) = —11-3 
as the dividing line. There are 9442 sources in the low 
Md population and 16 768 in the high M^ group. The 
GRAMS grid was intentionally designed to cover a wider 
range of stellar and circumstcllar shell properties than 
are spanned by actual AGB stars, and while the best 
fit models for these sources do have non-zero mass-loss 
rates, the rates are so low as to be consistent with zero, 
and we do not consider these low mass-loss rates to nec- 
essarily be physical. 

We highlight these low-M^ sources in Figure 
plot shows the [3.6] magnitude vs. the log (Variability 
Period). C-rich candidates are plotted in red, O-rich 
sources with \og(Md) > —11.3 are in blue, and the low- 
Md O-rich sources are highlighted in green. Visible in this 
figure are the 5 parallel AGB pulsation seque nces, labeled 
4 -1 a nd D from left to right (Figure 1 from | Riebe l et al. 
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This 



2010[ ). As discussed in |Vassiliadis fc Woo d (1993), stars 
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Table 8 

Sequence D Comparison 



Classification 




Sequences 1-4 


Sequence D 


C-rich 




2906 (26%) 


845 (19%) 


O-rich (all) 




7894 (73%) 


3505 (80%) 


O-rich (log(Af d ) > - 


-11.3) 


5405 (68%) 


2438 (69%) 


O-rich (log(M d ) < - 


-11.3) 


2496 (31%) 


1067 (30%) 


All Sources 




10 800 


4350 



Note. — Comparison of the population of Seq uence D with 
the radial pulsation P-L sequences 1—4 (see |Riebel et al.|201u} . 
The O-rich AGB stars in sequence D have the same ratio of 
high to low Mrf as the population of the other 4 sequences, 
consistent with sequence D population being drawn randomly 
from sequences 1—4. However, Sequence D is noticeably poor in 
sources classified as C-rich as compared to the other sequences. 

are expected to begin their evolution along the AGB near 
the bottom of sequence 4, and generally evolve to longer 
periods and to brighter luminosities, finishing their lives 
pulsating in the fundamental mode (sequence 1). The 
mechanism behind sequenc e D, the 'long secondary pe- 
riod,' is c urrently unknown (Nicholls et al. 2009; Wood & 
Nicholls||2009 |. We would expect the least-evolved AGB 
stars to exhibit the lowest rates of mass-loss, and we 
see that the low-M^ sources defined above cluster at the 
bottom of sequences 3 and 4, giving us confidence that 
these sources indeed represent the least-evolved sources 
in our sample, consistent with bare photospheres. In- 
deed, the low-M^ population is clustered almost exclu- 
sively on sequences 3 & 4, with almost none on sequences 
2 & 1. There is essentially no such thing as a non-mass- 
losing AGB star on the fundamental-mode sequence. We 
note, however, that the half of our sample with well- 
determined MACHO periods (that appear in Figure 17 1 



constitutes only 20% of the total of our entire sample. 
A full 80% of the total dust injection from evolved stars 
into the ISM comes from sources too red to have been 
detected by the MACHO survey, lack quality variability 
information, and therefore are not included in this figure. 
In a future paper, we will present initial constraints on 
the variability of these sources, using the warm Spitzer 
mission. 

There is also a population of low-M^ sources on Se- 
quence D. Table [H] summarizes the relative numbers of 
each type of source on sequence D as compared to the 
other sequences in toto. Sequence D has the same M^ 
population distribution as Sequences 1-4 (~l/3 low-M<j). 
Sequence D is more O-rich than the radial pulsation se- 
quences. The C- to O-rich ratio among the sequence D 
sources is only 24%, compared to 37% in the radially 
pulsating sequences. Bootstrap analysis shows that the 
difference in C- to O-rich ratio between sequence D and 
the other sequences is not consistent with sequence D 
being a randomly drawn subsample of the radially pul- 
sating sequences. 

4.4. The Total Dust Mass Return to the ISM 



Figure 18 shows the cumulative mass return to the ISM 
from AGB and RSG stars, as a function of M^ (left panel) 
and luminosity (right panel). Both panels show the run- 
ning total for both C- and O-rich sources separately, and 
for both combined. The O-rich grid extends to both 
higher luminosities (because it is intended to cover RSGs) 





Table 9 






Total M,j by population 




Population 


Total M d (xlO _6 M Q yr~ x ) 


Percent of total 


All Sources 


21.1 ±0.6 


100.0% 


C-rich AGBs 


13.64 ± 0.62 


64.6% 


O-rich AGBs 


5.5 ± 0.2 


26.0% 


RSGs 


2.0 ±0.1 


9.4% 


Extreme AGBs 


15.7 ± 0.6 


74.2% 



Note. — Total of Md broken down by classification. Column 3 lists 
the fraction of the total evolved star dust mass injection to the ISM 
each population contributes. Note that the category "Extreme AGBs" 
is a subset of O-rich AGB and C-rich AGB (most extremes arc C-rich). 



and higher values of M^ than the C-rich grid, but the C- 
rich stars contribute greatly to the total dust injection 
to the ISM, accounting for over half the dust produc- 
tion but only 21% of the sample by number. The total 
AGB+RSG dust mass return, (2.1 ± 0.06) x 10" 5 M Q 
yr is dominat ed by a very small fraction of stars. In 
agreement with Paper I and Boyer et al. (2011), we find 
that the total contribution or dust to the ISM from AGB 
stars is dominated by a very few 'extreme' stars. There 
are 1340 stars in our sample classified as 'extreme' using 
CMD cuts. Of these, 1299 (97%) are classified as C-rich 
by the GRAMS grid. These C-rich 'extreme' stars com- 
prise only 4% of our sample by number, but account for 
75% of the total Md- The few 'extreme' stars classified as 
O-rich are also heavy mass- losers, making up only 0.1% 
of our sample by number but representing 13% of the 
total dust-injection rate. We find a total AGB dust pro- 



duction rate of (1.91 ± 0.06) x 10~ 5 yr" 1 . This i s 



quite consistent with the work of Srinivasan et al.| ( |2009[ ) 
Table [9] breaks down our estimates of total M^ for each 
sub-population in our sample. We obtain exactly the 
same value for the du st injection rate due to RSGs as 



Matsuura et al | ( |2009[ ), 2 x 10~ b M o yr" 1 . Our results 
indicate that the dust production in the LMC is dom- 
inated b y C-rich AGBs, but not to the same extent as 



found by |Matsuura etTaC] ( |2009[ ) . 

Our estimate of the global dust production of the AGB 
p opulation in the LMC is a factor of ^2 smaller than that 



of Matsuura et al. ( 2009 1 . This discrepancy may in part 
be due to the systematically lower dust mass-loss rates 
in GRAMS carbon-star models than in the models used 
in that work 



(see § 5.1.2 



P aper V ) 

mill 



Additionally, re- 
indicates that SNe 



cent work by Matsuura et a 
may produce substantially more dust than previously es 
timated, up to ~1 M Q each. If this result is confirmed, 
our downgrade of the dust contribution of AGB stars 
serves to indicate that SNe may be the dominant source 
of dust in the ISM of the LMC. We would like to note, 
however, that care should be taken when comparing our 
results to other studies. For instance, the GRAMS grid is 
constructed with the assumption of a constant dust shell 
expansion velocity v cxp , which is kept fixed at 10 km s _1 . 
The actual value of t> CX p for LMC stars is uncertain, and 
it may depend on the luminosity , metallicity and gas: dust 
ratio (see, e.g. |van Loon||2000[ ) . The uncertain value of 
the expansion velocity translates to an uncertainty in the 
absolute value of the dust mass-loss rate and, therefore, 
the global dust injection rate. 
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Figure 18. Left Panel: Cumulative vs. M,; for individual sources. Dashed line: C-rich sources, dotted line: O-rich sources, solid 
line: total of all sources. The O-rich sources dominate the lower M,j region through their larger population. The C-rich sources are less 
numerous, but dominate at the largest values of M^. Right Panel: Cumulative vs. bolometric luminosity of individual sources. The 
same linestyles are used as in the left panel. Here O-rich sources, which contain both hot-bottom burning stars and red supergiants, extend 
to much higher luminosities, but it is the C-rich sources that dominate much more strongly. 



Table 10 

C-rich AGB M d vs. IR color 



Color 


Po 


Pi 


P 2 


Residual 


J - H 




-6.32 


0.56 


-6.01 


0.15 


H - [4.5] 


-24.50 


4.52 


-5.67 


0.09 


H - [c 


.01 


-28.18 


4.89 


-5.70 


0.08 


H - 241 


-29.43 


5.24 


-5.88 


0.09 


K B - 


3.6] 


-9.72 


1.91 


-6.10 


0.11 


K B - 


4.5] 


-13.94 


2.88 


-5.97 


0.12 


K B - 


5.8] 


-17.90 


3.64 


-5.88 


0.15 


K B - 


8.0] 


-15.40 


2.95 


-6.16 


0.09 


K B - 


24] 


-18.33 


3.59 


-6.25 


0.12 


[3.6] - 


[8.0] 


-9.37 


1.79 


-5.85 


0.15 



Note. — Hyperbolic functions of IR colors as 
observable proxies for in C-rich AGB stars. 

Listed are the coefficients of the best fit hyper- 
bola, of the form log M d - (col< ^ +P + P2. The 
'residual' column lists the MAD of the residuals 
to the best-fit hyperbola, a measure of the spread 
of the data about the best fit curve. 



4.5. Observational Proxies for Md 

Because of the observational investment necessary to 
obtain 12 bands of photometry for modeling pur- 
poses, or precise CO line profiles to trace the gas return 
to the ISM directly, it is desirable to develop simple, eas- 
ily observable proxies for this astrophysical quantity. We 
fit hyperbolas to Md as a function of all possible IR col- 
ors in our dataset. Not all colors serve as usef ul proxies 
for Md and we present selected fits in table 



10 



We have 

fit our functions to only the C-rich AGB population in 
our sample, as we find that the O-rich ACB and RSGs 
populations are very condensed in color space, and do 
not follow a well-defined empirical relationship. The pri- 
mary signature of circumstellar dust in O-rich AGB stars 
is the silicate feature at ~9 ^m, which is not adequately 
sampled by the IRAC camera aboard Spizter. The [9] /im 
band provided by the AKARI satellite is uniquely posi- 
tioned to probe this spectral feature, we hypothesize that 
K s — [9] would be an excellent proxy for in O-rich 
AGB stars, and intend to investigate this in future work. 
Figure 19 shows the relationship between and two 
[8.0] and [3.6] - [8.0]), derived from our 



colors (K t 



sample. We have chosen th ese colors for easy com pari- 
son to other studies, na mely Matsuura et al. ( 2009 ) and 
Gullieuszik et al. ( |2012[ ). In order to derive the combined 
dust budget for the LMC 
two-band color p roxies for Md developed by 



C, Matsuura et al. (20091 used 



gen et al. (2007). The left panel of Figure |19| compares 
the K s — [8.0] color proxy used in that work (dashed 
blue line) to the best fit relation to our data (solid green 



Groenewe- 



line). Our relation is detailed in Table 10 The right 
panel of Figure [19] shows another example of an IR color 
as a proxy for Md from C-rich AGB stars. Again the solid 



green line is the best fit to our data, detailed in Table [10 



Th e blue dashed curve is the proxy used by Matsuura et 
al. ( 2009[) , and the do tted blue line is the relation derived 
15y (jullieuszik et al.] ( |2012[ ) drawing on the early results 
of the IR VISTA surveyTjFthe LMC. Note that the curve 
plotted here has been divided by the ga s to dust ratio, 
ip — 200, used by Gullieuszik et al. (2012), who give their 
results in terms of total gas+dust mass return. This al- 
lows us to comp are the mass-loss r ates o f only dust. The 
method used by Gullieuszik et al. (2012) is similar to our 
own, and the agreement bet wee n our curve and theirs in 
encouraging (sec also Figure 13). The offset between the 
two curves is a result of different optical constants used 
to model the opacity of amorphous carbon dust used by 
our two teams. As our dust is more opaque, a given IR 
flux implies less total dust in our models compared to 
theirs, which corresponds to a lower M<j. There is cur- 
rently no clear consensus as to which set o f dust opacities 
is more correct. See § 5.1.2 in Paper V (and references 



therein) for a thorough discussion ot this effect 
5. CONCLUSIONS 

We present the largest sample of evolved stars sub- 
jected to radiative transfer modeling to date. By fitting 
^30 000 stars with ^12 bands of photometry each to pre- 
computcd radiative transfer models from the GRAMS 
model grid of O- and C-rich AGB stars and RSGs, 
we compute individual bolometric luminosities and dust 
mass-loss rates for the entire evolved star population of 
the LMC. Our work forms a useful prediction to be tested 
against models of stellar population synthesis and stellar 
evolution. 



We establish that the GRAMS model grid (Paper IV 
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0] color for C-rich sourc es. The solid curve in g reen is the best-fit line to the C-rich stars: 

12009) (see text for details). Right panel: Plot of 

9.37 
[8. ()]) + !. 79 ' 



Figure 19. Left panel: Plot of M<j vs. K B 

log Affj = ^ K ~^ o'])+2 95 — 6-16. The dashed line is the relation used by 

M^ vs. [3.6] — [8.0] color for C-rich sources. The solid curve in green is the best-fit line to the C-rich stars: log 



Matsuura et al. 



-5.85. 



The dashed curve is the relation used by |Matsuura et al.|^2009) , and the dotted curve is that derived by |GuIrieuszik et a.1. | ( [2012^ (see text 
for details). 

~2. That is, Carbon stars inject the same amount of 
mass into the ISM as O-rich AGBs, but two and a half 
times as much dust. 



Paper V ) is an important new tool for the study of stel- 
lar populations. Capable of generating statistically ac- 
curate predictions for difficult-to-observe stellar parame- 
ters such as mass-loss rate and luminosity for thousands 
of stars in a reasonable period of time, GRAMS will be 
valuable for int erpre ting the infrared observations of fu- 
ture missions (§ 2.3). Through comparison to previously 
published studies of AGB stars in the LMC, we show 
that the GRAMS O- and C-rich identifications are over 



80% accurate (§4.1) 



We present a quadratic formula for the K s band bolo- 
mctric corre ction of AGB stars as a function of J — K s 
color (§[42}. Our fit is detailed in Table [7j 

We characterize the luminosity functions of both O- 
rich and C-rich AGB stars in the LMC. The O-rich stars 
follow a distribution between Mboi = -3.5 - -5 centered 
at ~ -4.3 (~3500 L ). Individual C-rich AGB stars tend 
to be brighter, distributed between Mboi = -4 — -6 Mboi, 
centered at -4.7 Mboi (^5500 L ). Despite tending to be 
less luminous on an individual basis, O-rich sources domi- 
nate the integrated IR light from the LMC evolved stel lar 
population, by virtue of their greater numbers 



4.2) 



We find that the AGB stars exhibiting the 'long sec- 
ondary period' phenomenon, and lying on sequence D are 
significantly more O-rich than the stellar population on 
the other, radially pulsating sequences. We see a quali- 
tative relationship between M^ and pulsation period for 
radially pulsating AGB stars, in that sources with no de- 
tectable mass-loss are concentrated on sequenc es 3 & 4, 



and do not appear at all on sequences 1 & 2 (§ 4.3 ) 



We also derive hyperbolic fits to the M^ from C-rich 



AGB stars as a function of various near-IR colors (§ 4.5 ) 
Presented in Table [T0| our fits are derived from a sample 
much larger than those used by previous studies. 

We derive a total dust mass-injection rate for the entire 
RSG+AGB stellar population of (2. 13 ±0.02) x 10~ 5 M 
yr _1 of the LMC through consistent, direct summation. 
Assuming a gas-to-dust ratio ip = 500 for O-rich AGB 
stars and RSGs, and a ratio -0 = 200 for C-rich AGB 
stars, this translates to a total mass injection rate into 
the ISM from RSGs and AGB stars of M = 6.5 x 
10~ 3 M Q yr _1 , with half the mass coming from Carbon- 
rich AGB stars. Dominated by the uncertainty in the 
gas-to-dust ratio, this figure is accurate to a factor of 
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